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Understanding the mechanism of W-doping induced reduction of critical temperature ( T c ) for V0 2 
metal-insulator transition (MIT) is crucial for both fundamental study and technological application. Here, 
using synchrotron radiation X-ray absorption spectroscopy combined with first-principles calculations, we 
unveil the atomic structure evolutions of W dopant and its role in tailoring the T c of V0 2 MIT. We find that 
the local structure around W atom is intrinsically symmetric with a tetragonal-like structure, exhibiting a 
concentration-dependent evolution involving the initial distortion, further repulsion, and final stabilization 
due to the strong interaction between doped W atoms and V0 2 lattices across the MIT. These results directly 
give the experimental evidence that the symmetric W core drives the detwisting of the nearby asymmetric 
monoclinic V0 2 lattice to form rutile-like V0 2 nuclei, and the propagations of these W-encampassed nuclei 
through the matrix lower the thermal energy barrier for phase transition. 

Monoclinic/rutile vanadium dioxide V0 2 (Ml/R), a model metal-insulator transition (MIT) material, 
displays a first- order transition from a monoclinic insulating phase to a tetragonal metallic phase at a 
critical temperature (T c ) of 341 K, accompanied by a reversible switch in optical transmittance/reflec- 
tance in the infrared range 1 . This dramatic phase transition makes V0 2 potentially useful for 'smart' optical and 
electrical switches, which yet requires the regulation of T c to around room temperature (RT) 2 " 6 . The most 
effective route for appreciable reduction of T c is doping some donor impurity atoms, which could be easily 
realized by chemical substitution methods 7 " 10 . In particular, doping the V sublattice with W is very attractive since 
it yields a pronounced reduction in T c by —20 K/at.%W for the bulk and by —50 K/at.%W for the nanostructure, 
and even metallic ground state is realized for W concentration of >8 at.% 7,8 . Therefore, W-doped V0 2 system has 
been widely used as key building blocks for some smart devices, such as thermochromics and variable-reflectance 
mirrors 8 ' 9 . However, there are still many intriguing questions to be resolved: why such a small amount of W 
dopant could introduce so large a reduction in T c and what is the underlying physical mechanism to drive the 
MIT transition when W is involved? Right answers to these questions would enrich our understanding on the 
exact interplay of electron and lattice of correlated materials during MIT mechanism. 

Several experimental investigations have been carried out on the W x Vi- x 0 2 systems and attempted to present 
some structural and mechanistic explanations 10 " 12 . By using small-angle X-ray scattering and X-ray absorption 
methods, Tang et al. considered that the charger transfer from W ions weakens the V— V pairs, promoting the 
structural transformation from monoclinic to rutile 11 . Also, Booth et al found that a significant expansion in 
certain crystal axes induced by W dopants breaks the Peierls-paired V ions, leading to the reduction of T c 12 . 
However, these studies are limited to heavily-doped samples with relatively lower T c , in which the V0 2 host 
lattice is close to rutile phase at RT. As a result, only the final state affected by W doping is obtained, which can be 
found from their stable and almost invariable local structure (LS) around W atoms. On the other hand, the 
interatomic distances of W dopants are short in heavily- doped samples; therefore the strong coupling effect 
between the neighboring W atoms makes it difficult to isolate the structure and contribution of a single W atom 
from its environment. Therefore, the real W local structure and its effect on the surrounding monoclinic V0 2 
matrix, as well as its evolutions during MIT remain obscure. In order to attain physical insights into W-induced 
T c reduction, one promising way is to conduct detailed structural characterizations on specimen with wide 
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concentration gradients, which exhibits the reduction of T c from 
several to dozens of degrees. For this purpose a rather suitable tool 
is X-ray absorption fine structure (XAFS) spectroscopy thanks to its 
element selectivity and good sensitivity to local structures 13 , which 
allows one to separately probe and compare the LS around V and W 
atoms in a same sample, even in the dilute doping case. In fact, the 
capability of XAFS has been demonstrated in a recent work that 
provides a very distinct, qualitative signature for the monoclinic - 
to -tetragonal transformation in V0 2 14 . 

In this work, to completely understand the T c reduction mech- 
anism induced by W- doping, we study systematically the LS around 
W and V atoms in single crystals W x Vi- x 0 2 with a wide doping 
range from dilute to heavy concentrations (x= 0.54%, 0.95%, 1.45% 
and 2.38%) using the XAFS technique. These representative concen- 
trations allow us to present for the first time: (i) a clear figure of LS 
around an isolated W dopant in the monoclinic V0 2 lattice, and (ii) 
the reconstruction of concentration-dependent evolutions of W and 
V structures across the V0 2 MIT. 

Results 

The T c and the crystalline structures for the W x Vi- x 0 2 samples are 
determined by DSC and XRD as shown in Figure 1. DSC analysis 
indicates that increasing W content x from 0.54% to 2.38% signifi- 
cantly reduces T c from 54.1 to — 50°C (Figure la). The plot of T c as a 
function of W concentration in Figure lb shows an approximate 
linear relation with a slope of —20 K/at.%W, comparable with that 
~22 K/at. %W for thin films reported previously 7 . The RT XRD 
patterns displayed in Figure lc exhibit two characteristic diffraction 
peaks (asterisks) assignable to the monoclinic phase for the x=0.54% 
sample, suggesting the primary monoclinic V0 2 structure. At higher 
doping level (x> 1.45%), only the diffraction peaks corresponding to 
tetragonal phase can be observed 15 , demonstrating the complete 
transformation into rutile phase of V0 2 . It should be noted that 



the magnitude of the phase transition is weakened with increasing 
W concentration, reflected by the peaks in DSC results, the phenom- 
enon of which is well known in heavy atoms (W, or Mo) doped 
V0 2 716 . Although it is truly a practical problem, the W-doping is 
known as the most effective route to reduce T c , and has been con- 
sidered as the prototypical model in understanding the mechanism 
by which the T c of V0 2 MIT could be manipulated. 

To obtain detailed structural information, we performed element- 
specific XAFS measurements on both host V and doped W atoms. 
Figures 2(a) and 2(b) show the V K-edge extended XAFS (EXAFS) 
oscillation /(/c) functions and their Fourier transforms (FTs), 
respectively, along with the calculated EXAFS data for V0 2 with 
monoclinic (Ml, space group P21/c), rutile (R, space group P42/ 
mnm) and another monoclinic (M2, space group C2/m) phases as 
the references. With increasing W content, the LS around V atoms 
gradually evolve into that in tetragonal phase, as indicated by the 
featured changes in the k range of 4-8 A" 1 (Figure 2a) 14 . Correlated 
changes can be identified in FTs profiles (Figure 2b). The V0 2 lattice 
structure of 0.54% W doped sample remains in the monoclinic phase, 
because its FTs curve keeps the main features with that of pure 
monoclinic phase. When W content increases to 0.95%, the two V- 
O split peaks merge into a single peak at 1 .48 A and the V-Vl peak at 
2.12 A disappear, implying that the initial monoclinic V0 2 phase has 
been largely distorted towards tetragonal phase at RT. For the 1.45% 
and 2.38% W doped samples, the FTs curves match well with that of 
calculated tetragonal V0 2 phase. The structural transition from 
monoclinic to tetragonal phase with increasing W concentration 
clearly confirms the reduction of T c , consistent with the XRD and 
DSC results. 

Of note, the formation of the metastable M2 phase that has partial 
dimerization of the V-V chains was reported in the monoclinic to 
rutile phase transformation 1015,17 21 . However, metastable M2 phase 
is indeed tightly associated with the strains, and are nearly all found 
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Figure 1 | The T c and the overall structures for the prepared samples W x V!_ x 0 2 . (a) The DSC measurement results for samples with different W 
doping concentration, (b) The plot of transition temperature T c as function of W composition x for W % Vi_ x 0 2 single crystals. T c values are effectively 
reduced by W doping at —20 K/at.% W. (c) The XRD patterns measured at room temperature for W x Vi- x 0 2 samples. The diffraction peaks correspond 
to the rutile phase, while the additional peaks (indicated by asterisks) correspond to the monoclinic phase. 
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Figure 2 | XAFS results for W x V!_ x 0 2 . (a) The V iC-edge EXAFS oscillations and (b) their Fourier transforms for W x V!_ x 0 2 samples, along with the 
V0 2 theoretical results of monoclinic (Ml), rutile (R) and another monoclinic (M2) phases for references, (c) The W L 3 -edge EXAFS oscillations and (d) 
their Fourier transforms, along with the reference W0 3 result, and the simulated curves for W-substituted monoclinic and tetragonal V0 2 phases. Fitting 
results are plotted as open circles. 



within the nanostructures that have strong surface-induced stress to 
stabilize M2 phase 10,1517 , while there is no evidence for isolation of 
stable M2 phase for large V0 2 crystallites without apparent stresses 
and defects 18 . Even in doped bulk V0 2 crystallites, only the Al and Cr 
doping that could result in the increase of T c are reported to stabilize 
the M2 phase 19 20 . Compared with W ions, these dopants have smaller 
ionic radii and are acceptor impurities, which would lead to the 
differentiation of two types of V atoms, corresponding to the pro- 
posed M2 monoclinic phase 21 . Hence, we consider that the inter- 
mediate M2 phase can be precluded in our highly crystalline W- 
doped V0 2 , in agreement with the previous reports 12,22 . This conclu- 
sion can be further confirmed by the significant differences in the 
oscillation x(k) and FT curve between the calculated spectrum of M2 
structure and our samples. 

We further focus on the LS around W atoms that could be deduced 
from W L 3 -edge XAFS data (Figures 2c and 2d). The LS around W 
atom should be symmetric in nature similar to the reference W0 3 . 
This is inferred from the fact that these samples, like W0 3 , also 
present two peaks at 4.0 and 6.2 A" 1 in the low k region (3-8 A" 1 ) 
of the k%(k) oscillation curves (Figure 2c), and only one prominent 
W-O peak without splitting appears in the FTs curves (Figure 2d). It 
should also be noted that the symmetry around W core at the low 
concentration (< 1%) is distorted due to the significant compression 
of almost isolated W dopants by the overwhelmingly asymmetric 
monoclinic V0 2 lattices, which can be reflected by the missed fine 
features (marked by arrows) in the 3-8 A" 1 region, as well as by the 
damped oscillations in the 8-14 A" 1 region. Furthermore, we have 
also calculated the W L 3 -edge FTs spectra for the W-substituted 
monoclinic and tetragonal V0 2 phases. Interestingly, the FT curve 
of low-doped (0.54%) sample is more similar to the calculated 
spectrum of the W-substituted tetragonal phase than that of 



W-substituted monoclinic phase which shows an obvious splitting 
for W-O peaks. Only a slight difference can be found for the position 
of W-O peaks which shifts from 1.35 A of W0 3 toward 1.51 A of W- 
substituted tetragonal phase with increasing W concentration. 
Although the host V0 2 is in monoclinic structure inferred by XRD 
results, the local environment of W dopant in the monoclinic V0 2 
matrix has yet assembled the tetragonal-like structure with com- 
pressed W-O pairs, implying the W-doping induced structural 
transition of the neighboring host V0 2 lattice. Increasing x, the 
W-Vl and W-V2 peak positions (2.42 and 2.97 A, respectively) 
increase to 2.67 A and 3.13 A, respectively, along with the enhance- 
ment of their intensities, implying the expansion of W local struc- 
tures. The FT curves for highly doped samples (x>0.95%), 
particularly the W-V peaks, are almost independent of the W con- 
centration, and this trend continues until the heavy doping level of 
5% < x <33% as reported in Ref. 11. These indicate that the LS 
around W atoms would not undergo remarkable changes when the 
V0 2 host lattice is already in or near the rutile phase if the W con- 
centration reaches a certain level, confirming the intrinsic and stable 
nature of symmetric W structure. 

The quantitative structural information was extracted from 
EXAFS data fitting using the ARTEMIS module of IFEFFIT. The 
extracted parameters are summarized in Figure 3 and Table I and 
II. The obtained V-O and V-V distances show that the monoclinic 
V0 2 host lattice is largely distorted in the x=0.95% sample, and has 
completely transformed into tetragonal rutile structure when x 
increases to 1.45%. Because of the slight distortion of W atoms by 
the surrounding V0 2 lattice, the EXAFS data for x=0.54% sample 
were fitted using the monoclinic phase with degenerate W-O bonds. 
The rest three can be fitted well by using tetragonal rutile phase, due 
to their higher symmetric structures. The trend in Figure 2b shows 
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Figure 3 | The bond distances as a function of W concentration. The 

bond distances of the nearest three shells of (a) V and (b) W as the center 
atom, based on the XAFS fitting results of each sample with different W 
doping concentration. 



steady increase of the average W-O and W-V bond distances, con- 
firming the gradual expansion and relaxation of LS around W atoms 
with increasing W content. Here, we have found the expansion of 
both W-O and W-Vl pairs, which however are shown constant in 
Ref. 12, with the W concentration. The whole expansion observed 
here can be ascribed to the distinct W concentration gradient, thus 
reflecting the real evolution of LS of W dopants in monoclinic V0 2 
across MIT. 

Discussion 

Summarizing our experimental results, we schematically show in 
Figures 4a~c the crystal structure variations for W-doped V0 2 
samples with three representative concentrations (0.54%, 0.95% 
and 2.38%). For the purpose of clarity, we magnify the structural 
model around W atoms in a planar visual. At first sight, the LS 
around W atom highlighted in white exhibits two prominent fea- 
tures: the increased symmetry and the expanded volume. A point 
worth emphasizing is that several XAFS studies by Patridge et al. 23 , 
Booth et al 12 , and Whittaker et al. 10 , as well as the static lattice 
calculations by Netsianda et al. 24 , had suggested that dopants induce 
a local structural distortion towards high symmetry of rutile phase, 
which is generally accepted in the work of doping-induced phase 
transition of V0 2 . Moreover, our results not only confirm this basic 
model, but also demonstrate a concentration -dependent evolution of 
LS of W atoms. We find that the LS of W atom at low concentration is 
tetragonal-like but distorted by surrounding asymmetric monoclinic 
V0 2 lattice; at high doping level, the repulsion of these distorted 



tetragonal units drives them to restore and stabilize the symmetric 
structure. This structural evolution of W atoms reflects the strong 
interaction between W and V0 2 structures that triggers the struc- 
tural derivations of near V0 2 lattice from asymmetric monoclinic 
phase towards more symmetric. 

To reveal the reason why a small amount of W doping lowers the 
T c , we performed density functional theory (DFT) calculations with 
Vienna Ab-initio Simulation Package (VASP) code based on the 
atomic structural parameters by XAFS. The calculations employed 
the projector augmented wave (PAW) method and the generalized 
gradient approximation (PBE-GGA) for electron exchange-correla- 
tion interaction 25,26 . A super large 4X3X5 supercell with 432 atoms 
containing 143 V atoms, 288 O atoms and 1 W atom, correspond- 
ing to —0.6 % W concentration was used. At this concentration, the 
V0 2 lattice is of the monoclinic structure in the temperature range 
from 0 K to RT. We first considered the spatial extension of the 
influence on the atomic and electronic structures of V0 2 by W dop- 
ing. Figure 4d shows the relaxed structure projected along [010], in 
which the distortion degree of V0 2 lattice is characterized by the 
difference of bond length \ AR\ (\AR\ = liWorted - Venous I )• We can 
find that the structural change shows a dopant- centred nature, since 
I AR| near W dopants are all larger than 0.01 A whereas those far away 
from W atoms are <0.002 A. Furthermore, we examine the spatial 
redistribution of charge density upon doping W ions (Figure 4e). 
Also, the electron density pile-up is only significant in the neighbor- 
hood of W atoms along or cross the c axis. This is reasonable because 
the V0 2 structures away from W atoms remain in the insulating 
monoclinic phase which inhibits electron mobility 27 . Importantly, 
the relaxed structures in the theoretic calculations confirm that W 
dopants could induce large variations in both atomic and electronic 
structures of nearby V0 2 lattice. In this sense, the corporation of local 
lattice disorder induced by the insert of W ions with larger volume, 
and the pile-up of itinerant electrons on the neighboring V-sites 
released by W ions with higher valence suppress the dimerized lattice 
distortion and destabilize insulating phase of V0 2 . 

On the basis of experimental and theoretical results, we consider 
that the W-doping effect is confined in a local area (covered by each 
blue area of W cores schematically illustrated in Figure 4f), in which 
the V0 2 has accomplished the phase transition, and can be seen as 
the nucleus for MIT. Due to the existences of these pre-MIT nuclei, 
the thermal energy needed to trigger MIT is lowered compared with 
that for pure V0 2 single crystals 28 ' 29 . At the low doping concentra- 
tion, the nucleus can be seen as isolated, while at higher concentra- 
tion two neighboring nuclei could touch each other due to the 
increased average W-W distance (the overlapped area shown in 
Figure 4f). Meanwhile, the XAFS results of the steady expansion of 
W structures and significant phase transition of V0 2 lattices sug- 
gested the propagation of W-induced changes towards the farther V 
atoms. Therefore, the increased formation of nuclei, along with their 
percolation in the materials, will further lowers thermal energy to 
fulfill the whole phase transition as expected. This picture of concen- 
tration scaled influence of W dopants is somewhat like the "metallic 
nanopuddles" in thermal-induced pure V0 2 MIT process 30 , in which 
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Qazilbash et al. have revealed that nanoscale percolative MIT took 
place first with metallic puddles in insulating host, and subsequent 
heating drives its complete transformation to metallic phase. 

In conclusion, the behavior of an isolated W dopant and its influ- 
ence on the V0 2 host lattice are for the first time figured out by X-ray 
absorption spectroscopy coupled with the first-principle calcula- 
tions, which helps us to unravel the mechanism of the reduction of 
MIT T c induced by W-doping. A clear atomic structure of isolated W 
dopant has been unveiled, which possesses an intrinsically symmet- 
ric tetragonal-like structure, driving the detwisting of the nearby 
asymmetric monoclinic V0 2 lattice towards rutile phase. These W- 
encompassed rutile -like V0 2 regions in fact act as nuclei for the 
"metallic puddles" which propagate through the host matrix, thus 
lowering the thermal activation energy across phase transition. These 
results also provide a clear indication to modulate the T c by doping 
impurity with symmetry corresponding to that of targeted structure, 
which has been noticed in some syntheses. 



Methods 

Synthesis and Characterization. Single crystals of W x Vi_ x 0 2 used in this study were 
prepared by reacting vanadyl acetylacetonate (VO(acac) 2 ) with the required amount 
of sodium tungstate (Na 2 W0 4 ) in aqueous solution at 200°C for 12 hours. The 
resultant products were calcined under argon at 700 °C for 8 hours. The phase 
homogeneity and the reversible transition characteristics of the samples were 
confirmed by X-ray diffraction (XRD) and differential scanning calorimetry (DSC). 
XAFS data at the W L 3 and V K-edges were collected in fluorescence and transmission 
modes, respectively, at ambient temperature (25°C) on beam lines U7B and U7C of 
the National Synchrotron Radiation Laboratory, China. The storage rings of NSRL 
were operated at 0.8 GeV with the maximum current of 200 mA. The hard x-ray 
beam was from a three-pole superconducting Wiggler with a magnetic field intensity 
of 6 T. The fixed-exit Si(lll) flat double crystals were used as monochromator. 
Ionization chambers filled with Ar/N 2 , and a 7-element pure Ge solid-state detector 
were used for the transmission and fluorescence x-ray detection, respectively. 

XAFS data analysis. The EXAFS signals, %(k), were extracted from the absorption 
raw data, |i(E) from the averaged experimental spectra of four repeated scans, by 
using the ATHENA module 31 implemented in the IFEFFIT package 32 and the 
USTCXAFS software packages 33 . The %(k) data were weighted by k 3 and k 2 for V and 




Figure 4 | Models and theoretical calculations. Models of LS distortions of three representative doped concentrations (a) 0.54%, (b) 0.95%, and (c) 
2.38%. The white area displays the LS evolutions around W atoms. The variations of average W-W distances and [W0 6 ] octahedral geometries are shown 
in the left top and bottom panel, respectively, (d) The relaxed structure projected along [010] axis and (e) the corresponding spatial redistribution of 
charge density of W doped V0 2 with about 0.6% concentration; I AR\ ( I AR\ = I ^distorted - ^previous' ) based on calculation results give a distinct pattern of 
doped effects. The blue and red isosurfaces correspond to charge densities of 0.005 e/A~ 3 and —0.005 e/A" 3 , respectively, (f) A schematic diagram of 
concentration scaled influences of W dopants in V0 2 , in which shows the low concentration with isolated dopants and the high concentration with 
interactional dopants. 
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W, respectively. The normalized absorption data were Fourier transformed using a k- 
space range of 2.8-14 A" 1 . The theoretical scattering amplitude and phase shift 
functions for all the paths were calculated by using the FEFF7 code. The V K-edge 
XAFS spectra below and above 0.95% W concentration can be fitted well using the 
structural model of monoclinic V0 2 (P21/c) and tetragonal V0 2 (P42/mnm) based 
on ICSD #15889 and #4110, respectively. According to the T c , W L 3 -edge XAFS 
spectra below and above 0.95% concentration was fitted in the framework of 
monoclinic and tetragonal phase, respectively. The FEFF input files for the fitting of 
W L 3 -edge XAFS spectra was altered by replacing the V core in the monoclinic and 
tetragonal FEFF input files with W. The fitting was done in the R range [0.8, 3.8] A for 
all the data. In order to decrease the number of independent parameters, in the fitting 
procedure the coordination number was fixed to the nominal value for each scattering 
path, because the coordination numbers are not changed during MIT. For the 
monoclinic phase, the variation of M-O (M= V and W) bond distances in the range of 
1.76-2.03 A are classified into three groups, i.e., (1.76 A), (1.86 A, 1.87 A), 
(2.00 A, 2.01 A and 2.04 A) according to their contribution to the scattering paths, and 
four kind M-V2 bond distances (3.36 A, 3.46 A, 3.56 A and 3.66 A) are included in the 
fitting. For each shell, the interatomic distance R, the Debye- Waller factor a 2 , and the 
edge-energy shift AE 0 were allowed to vary. 
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